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Abstract 



O I We study a mean-field model of a Kondo alloy using numerical techniques 

and analytic approximations. In this model, randomly distributed magnetic 

(N' 

^ ' impurities interact with a band of conduction electrons and have a residual 

o\ • 

en ' 

RKKY coupling of strength J. This system has a quantum critical point at 

(N 

O ' J = Jc ^ T^, the Kondo scale of the problem. The T dependence of the spin 

G> ' 



susceptibility near the quantum critical point is singular with x(0) ~ xi^) oc 
T^ and non- integer 7. At J^ 7 = 3/4. For J ^ Jc there are two crossovers 



'^ ' with decreasing T, first to 7 = 3/2 and then to 7 = 2, the Fermi-liquid 



value. The dissipative part of the time-dependent susceptibility x"{^) oc a; as 
uj ^ except at the quantum critical point where we find x"{^) ^ V^- The 

X" 

Vh ' characteristic spin-fluctuation energy vanishes at the quantum critical point 

with ujsi ~ (1 — J/Jc) for J < Jc, and cOsi oc T^'^ at the critical coupling. 
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I. INTRODUCTION 

The understanding of metallic states which do not fit within the framework of Fermi- 
liquid theory is one of the important current challenges of condensed matter physics [|l|. 
This issue is relevant to a large class of /-electron materials which present anomalies in their 
thermodynamic and transport properties at low-temperature 0. Two important features 
characterize the physics of these systems. One, is the interaction of the conduction electrons 
with localized magnetic moments via the Kondo coupling. The other, is the inevitable 
presence of disorder due to the alloying process. Several models in which non-Fermi-liquid 
(NFL) behavior arises as a consequence of the interplay between disorder and magnetic 
interactions have been proposed in the literature. In the Kondo-disorder model of Miranda et 
al. 0, randomness in the local hybridization matrix element between localized and itinerant 
electrons is thought to be at the origin of NFL behavior. In this theory, the disorder generates 
a broad distribution of Kondo temperatures whose tail extends down to T^ = 0. Therefore, 
a finite fraction of the localized spins remain unquenched at all temperatures and gives rise 
to singularities in the thermodynamic and transport properties. 

In the metallic spin-glass model [§|-|^ , the focus is on the consequences of randomness in 
the RKKY intersite couplings. This type of disorder is modeled by including a spin-glass- 
like exchange term in the Hamiltonian. The system has a quantum phase transition when 
the strength of the magnetic interaction J becomes comparable to the Kondo temperature 
of the underlying Kondo lattice, T^ . Beyond this point, the ground-state is no longer 
a non-magnetic metal but it exhibits long-range spin-glass order. NFL behavior results 
from the power-law behavior found in the neighborhood of the quantum critical point, a 
scenario that is similar to that proposed to explain NFL behavior in systems close to ferro- 
or antiferromagnetic instabilities [0. 

In this paper, we study a mean-field Kondo-alloy model recently proposed and discussed 
by Sengupta and Georges |Q. In their paper, these authors did not solve the original 
Hamiltonian but a simpler solvable quantum-rotor problem [R] that was assumed to exhibit 



the same low-frequency behavior. Here, we solve numerically the Kondo-alloy model using 
classicaland quantum Monte Carlo techniques. We find a quantum phase transition at J = 
Jc K. 1.15 T]^ where the zero-temperature spin-glass susceptibility of the system diverges. 
At the critical coupling, the T-dependence of the uniform magnetic susceptibility is singular 
with x(0) — x(T) ~ T^/^. This anomalous T-dependence is also found above the crossover 
line T / Jc ~ (1 — J j Jf)'^^'^ . Below this line, we still find unconventional behavior but the 
exponent is different, x(0) — x(T) ~ 2"3/2^ p^^^ j _^ j^^ ^j^^ normal behavior, bx ~ T^, is 
recovered at low enough temperature. The numerical results for the frequency dependence 
of the susceptibility are very well described over a wide range of temperature and frequencies 
by a simple approximate expression that we derive from the original model in the strong- 
coupling limit. The strong coupling solution reduces to that of the simplified model of 
Sengupta and Georges in the u; — i> limit. The spin-fiuctuation spectrum is Fermi- 
liquid-like for cj — * everywhere except at J = Jc- We find x"(^) oc a; for cj ^ cUsf where 
the spin-fluctuation frequency cjsf oc (1 — J/Jc) for J ^ Jc and tUgf oc T^/^ at the critical 
coupling. At the quantum critical point, 'x"{uj) oc ^/uj. This implies a slow decay of the 
spin-spin correlation function, < Sz{t)Sz{^) >~ t^'^/^, that anticipates the appearance of 
long-range order in the system. 

The paper is organized as follows: In Section |l| we introduce the model Hamiltonian 
and use it to derive an effective local action for the spin degrees of freedom. In Section 



mi , we discuss two equivalent fromulations of the effective model that are well suited for a 
numerical investigation of the problem. These are based on the formal equivalence between 
the Kondo-alloy problem and two other models. The first one, is a classical one-dimensional 
Ising chain with short- and long-range ferromagnetic interactions, and may be solved by 
classical Monte Carlo simulation. The second model system describes a single S = 1/2 
quantum spin evolving in the presence of two magnetic fields, one that is time-dependent and 
random in the longitudinal direction, and another that is static and fixed in the transverse 
direction. This problem may be solved using a quantum Monte Carlo algorithm. The results 



of the simulations are presented in Section |^ where we also derive a simple analytical 



approximation which allows for a transparent interpretation of the data. This is followed by 
a comparison of our results to those obtained by other authors. 

II. THE MODEL 

In a disordered Kondo alloy, randomly distributed spins interact with a band of conduc- 
tion electrons through a local Kondo coupling. There is also a residual RKKY exchange 
interaction between the spins which is random because of the disorder in their positions. 
Many of the systems studied experimentally exhibit uniaxial anisotropy as a result of strong 
crystal-field and spin-orbit effects. Therefore, to a first approximation, only the coupling 
between the components of the localized spins along the easy-axis needs to be considered. 
The simplest model with these characteristics is a Kondo-lattice model with an additional 
Ising-like random exchange term. The Hamiltonian of the model is : 

H = -Y: ^clc,^ + J|f E St St + ^ E {S^sr + h.c.) -W J^^StS^. (1) 

i,j,cr i i i,j 

Here, Si is a localized spin operator at the i-th site of a lattice of size A^. The creation and 
destruction operators for the conduction electrons are c^ and Cjo- and Si = 1/2 J2a,i3 cla^a^isCifs 
is the local electronic spin density. The nearest-neighbor electron hopping integral is tij = 
tj \fz where z is the connectivity of the lattice, and J^ and Jj^ are the longitudinal and 
transverse Kondo couplings, respectively. The nearest-neighbor couplings between the spins, 
Jjj, are quenched random variables for which we assume for simplicity a Gaussian probability 
distribution with zero mean and variance (Jfj) = J'^/z. The z-dependent normalization of 
the off-diagonal couplings t^j and Jij has been chosen such that the results in the 2; — >■ cxo 
limit to be considered below are finite 0. 

The last term on the righthand side of Eq. |l] is the well-known Sherrington Kirkpatrick 
model [|10 that has a phase transition to a spin glass state at T° oc J. The local Kondo 



coupling favors screening of the localized spins by the conduction electrons below a charac- 
teristic temperature T^. As a consequence of the competition between these two terms, a 



spin-glass ground-state is only possible for J > Jc ^ T^- Therefore, Tg ^ as J -^ Jc from 
above and the system remains paramagnetic down to zero temperature for J < J^.. The 
point T = 0, J = Jc where the nature of the ground state of the system changes defines the 
quantum critical point |[7|JTl[]. 



We investigate the properties of model (|I|) near the quantum critical point in the frame- 
work of a dynamical mean- field theory |jl2| . In this approach, exact in the limit of infinite lat- 



tice connectivity, the degrees of freedom on any particular lattice site are isolated and treated 
exactly, while the rest of the system is replaced by an effective medium to which the chosen 
site is coupled. The properties of the effective medium are determined self-consistently from 
the solution of the single-site problem. In the limit z ^ oo, the configurational average over 
the random couplings can be performed explicitly [|1^] and the intersite terms in Eq. (|l]) can 



be eliminated. This reduces the problem to a magnetic impurity embedded in an electronic 



bath and subject to a dynamic magnetic self-interaction [J12|,|T3[. Ignoring for the moment 
the anisotropy of the Kondo coupling in order to simplify the notation, the effective action 
of the single site problem in the paramagnetic phase may be written as ^ : 

S\oc= - I / dTdT'cl{T)gQ^{T -t')c^{t') + Jk I dTS{r)-s{T) 
Jo Jo Jo 

J2 r/? rP 



drdr'S,{r)x{r-r')S,{T'), (2) 

z Jo Jo 

The functions xi'^) ^i^nd ^o(t), a priori unknown, are determined by the feedback effects 
of the coupling of the impurity site to rest of the system through a set of self-consistency 
conditions. Their precise form depends upon the shape of the non-interacting electronic 
density of states J\f{e) of the lattice [^ . In the case of a semicircular density of states, the 



self-consistency equations acquire a particularly simple form. We have : 

x(r) = (T(S.(r)S,(0)))5,., (3) 

for the magnetic degrees of freedom, and 

So\r - r') = (-^ + /.] 6{r - r') - fg{r - r'). (4) 



for the electronic degrees of freedom. In the above equations, T is the time-ordering operator 
along the imaginary-time axis < r < /3, /i is the chemical potential and 

^(r) = -(r(c(r)c+(0)))5,, (5) 

It follows that Q and x ^i^re, respectively, the exact local electronic Green function and 
the imaginary-time dependent spin susceptibility. For general A/'(e), Eq. (^ is replaced by 
a more complicated implicit condition [Q. 

The solution of this set of coupled self-consistent equations is still a very difficult task. 
It may be argued, however, that knowledge of the exact bath Green function is not essen- 
tial for the understanding of the low-frequency spin dynamics of the model. This follows 
from a perturbative argument ||^ that establishes that the long-time behavior of the exact 
bath Green function is qualitatively the same as that of a bath of non-interacting electrons, 
i.e., Qq{t) ~ 1/r. But the form of the low-energy effective action for the localized spins is 
determined precisely by the asymptotic behavior of the electronic Green function. There- 
fore, if we ignore Eq. (^) and fix Qoiiujn) = X!^ deAfo{e)/ {iun + fi — e) where A/o(e) is the 
unrenormalized density of states, we will still get a qualitatively correct description of the 
low-frequency spin dynamics of the model. 

Further progress can be made by performing a Hubbard- Stratonovich transformation 

that decouples the last term in Eq. (0) . Introducing a set of time-dependent random fields 

?7(r) that couple to the spin operators ||l5l , the partition function of the problem may be 

rewritten as 

1 rP r(^ 



Vriir) exp 



dTdT'rj{T)x ^{t -T')rj{T') 
2 Jo Jo 



Zk [v] , (6) 



where 



and 



ZkIv] = fVc{T)Vc+{T)Tr{rexp-SK), (7) 

J ^z 



Sk=- I I dTdr'ct{r)g^\r - r')c„{r') + Jk f drS{T) ■ s{t) 
Jo Jo Jo 

-J f dTS,{T)r]{T). (8) 



Eq. (§) is the action of a single Kondo impurity in a time-dependent magnetic field Jti{t) 
in the z-direction. Within dynamic mean-field theory the partition function of the Kondo 
alloy is thus given by the average over all the realizations of the random field of the partition 
function of the modified Kondo problem of Eq. (^). 

Equations (|§)-(|§) subject to condition (0) define the mean-field model of the Kondo alloy. 
In the next Section we shall show that this model may be cast in two different forms both 
of which are well suited for setting up schemes for the numerical solution of the problem. 

III. METHOD 

A. Formulation of the problem 

As we are only interested in the spin dynamics of the system, the electronic degrees of 
freedom in Eq. (^) may be integrated out. In the case of the single-impurity Kondo model 
this leads to the well known Coulomb gas representation [|r^ of the partition function of the 
problem. The same is true for the generalized problem of Eq. (^) as the additional random 
term commutes with the longitudinal part of the Kondo coupling. The Anderson- Yuval 
technique [|l^] may therefore be applied to Eq. (0). After averaging the resulting expression 



over the distribution of random fields, Eq. @ may be rewritten as : 

^CG = E/ dn dT2... dT,J^] exp ^(_iy+^F(r,-r,) , (9) 

where 

_Li = 2„^_j„„-._ + j.,(,). (10) 

and To is a short time cutoff of the order of the inverse bandwidth of the electron bath. 
The coupling constant a = (1 + ^)^ where 5 = — tan^^ (7rroJ|^/4j is the phase shift for 
scattering of electrons from a local potential J\f /4:. 

In Eq. (^, Tj, i = 1, . . . ,2n are the positions on the time axis of successive spin-flips 
generated by the transverse part of the Kondo coupling and the function V{t) represents 



the interaction between pairs of spin-flips. The first term on the righthand side of Eq. ([lOD 



is famihar from the work on the Kondo model [|T^]. It arises from the singular response 
of the conduction electron bath to a spin flip on the impurity site. The second term is 
characteristic of the alloy model and represents the reaction of the rest of the spins to the 
local perturbation. It is interesting to notice that the partition function of a recently studied 
extended two-band Hubbard model [l^ can be cast in a form equivalent to Eq. (|^). 

Equation (P) is not yet in a form suitable for computation of the magnetic correlation 
function. We shall next establish a formal equivalence between Zqg and the average partition 
function of a single quantum S = 1/2 spin in the presence of a random Gaussian time- 
dependent longitudinal field ^(r) and a static transverse field F, a problem that can be 
solved numerically using the quantum Monte Carlo method of references [jl5[ and ||18| . The 
partition function of the quantum spin problem Zqs is : 



^QS 



P.^exp 



X Tr T exp 



f fdrdr'i{r)Q-\r,r')i{r') 
Jo Jo 



1 f/^ ff^ 
2 





(11) 



where Q{t) is the correlation function of the random component of the magnetic field. 
To demonstrate the equivalence of @ and ([Til), we first perform a Trotter decomposition 



of the time-ordered exponential in Eq. (|Tl 

r/3 



Texp 



dTh{T) ■ S{t) 



M 



k=l 



Yl exp Ar/i(rfc) ■ S 



(12) 



where At = (3/M. We next introduce a complete set of intermediate states |o"fc)(o"fc| at each 
imaginary time-slice Tk- The matrix elements in the Trotter expansion may be evaluated 
using the expression 



{a\ exp {Ar Ut)S,{t) + TS,{r)]} \a') ^ e^^«(^)- 



FAt 



(13) 



valid in the limit Ar -^ 0. After averaging over the field ^(r) and taking the limit M — > cxd, 
the partition function of the model may be expressed as a sum of contributions of individual 
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°° p ll r0 ri^ 

^QS = E / ^^^"^ exp - / / dTdT'a^''\r)Q{r - r')a(")(r') + nlnT 



(14) 



"spin histories", each of these being one of the possible sequences of the eigenvalues cr(r) 
±1/2 of the intermediate states. We find : 

"1 /•/3 /-/S 
/ Va"-"' exp 

71=0 

where a^"-* (r) is a spin history with n spin flips in the interval < r < /5 and the integration 
is over their positions. Eq. (P) follows from Eq. (Q) by integrating twice by parts the first 
term in the exponential, provided that we choose T = J^ and that we identify Q{t) with the 
righthand side of Eq. (|lOD- The original problem has thus been reduced to the evaluation of 
the partition function of Eq. ( pT]) subject to the condition @. 

An alternative numerical method may be formulated by taking advantage of the asymp- 



totic equivalence |1T9| , |20[| between the Coulomb gas representation (g) of the Kondo prob- 
lem and the partition function of a classical one-dimensional Ising spin chain with nearest- 
neighbor and long-range interactions. This problem may be solved numerically using stan- 
dard classical Monte Carlo techniques as has been recently done for the single-impurity 
anisotropic Kondo model in reference [pO]. It may be shown by a straightforward general- 



ization of the methods of reference [jT9[ that the Ising-chain model relevant for our problem 



IS 



Zi = Y^ exp 

{s.} 



E Kt^nSiSi+i + E ^LR{i - j)SiSj 

i<L i<j<L 



(15) 



where the S'i = ±1 are Ising variables and the number of sites in the chain is L = /^/tq. 
The spin-spin interaction consists of a short-range part, K^n ~ — l/21n (Ji4'ro/2), and a 
long-range part given by 



K 



LR[i-J, 



2a {n/NY 



sin^ [7r(j - 1)/N] 



+ rr^xirolt - j\) 



(16) 



It is worth noticing that, while both of these approaches can be used to solve the present 
strongly anisotropic Kondo-alloy model, only the first one can be generalized to the case of 
a non-Ising spin-spin interaction. 



B. Numerical methods 

We have simulated the mean-field Kondo-alloy model using the two formulations de- 
scribed in the previous section as each has its own advantages and drawbacks. In particular, 
while the systematic error introduced by the discretization of the imaginary time is larger 
for the quantum simulations, statistical fiuctuations are far more important in the classical 
case. We have empirically found that the latter method is more accurate for the compu- 
tation of the static susceptibility at low temperatures whereas the former one gives better 
results for the overall frequency dependence of the spin correlation function. 

The numerical procedure used to solve the self-consistent problem is as follows : i) an 
approximation to x{^) is used as input in either Eqs. ([10|) or (pJBI) . ii) the spin-spin correlation 
function is obtained from a Monte Carlo simulation (see below), iii) a new xi'^) is computed 
from condition @ and used as a new input in step i). This procedure is iterated until 
successive values of the correlation function differ by less than a fixed tolerance level (see 
below). This takes from four to fifteen iterations depending on the temperature and the 
values of the parameters. 

The simulations of the classical problem defined by Eq. (pISj) have been done using a 
standard Monte Carlo heat-bath algorithm for Ising chains of up to 256 sites. The quantum 
problem of Eq. ([IT]) has been simulated as follows. The imaginary-time axis is discretized in 
slices of width At = j3/M and the time-ordered exponential appearing in Eq. ([llD is approx- 
imated by a Trotter product of M factors. The statistical weight of a given configuration 
{^(r)} is thus expressed in terms of the trace of a product of 2 x 2 random matrices. The 
corresponding contribution to the spin correlation function x(r, r') is computed by insert- 
ing two additional az Pauli matrices at the appropriate places in the matrix product. It is 
important to choose the parameter M appropriately. If M is too small, the systematic error 
introduced by the Trotter approximation is large. If M is too large, however, the algorithm 
is prone to numerical instability. We found that the choice j3 ^ 0.25Mro with M < 128 is a 
satisfactory compromise. This sets a lower limit to the temperatures that we can simulate. 
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To^min ~ 0.03. The simulation is most conveniently done in the space of the Masubara- 
frequency components of the field, ^{un) = Jq dT^{T) exp{—iuJnT) |]15|. These are finite in 
number as a consequence of the discretization of time : Un = 27mT with n = 0, . . . ,M — 1. 
In an elemental Monte Carlo move, a change of the complex field C,{uJn) for a single fre- 
quency is attempted. A full Monte Carlo step is complete when elemental changes have 
been attempted for all the Matsubara frequencies. 

The precision of the numerical calculations presented below is determined by two factors, 
namely, the statistical error of the Monte Carlo calculation and the stopping criterion used 
in the enforcement the self-consistency condition. A typical quantum Monte Carlo run 
consisted of 4 x 10^ Monte Carlo steps per time slice. This corresponds to an absolute error 
of the order of 2 x 10~^ in x(t). As we mentioned above, the simulations based on Eq. (p!5D 
are noisier than those based on Eq. ([TTD which requires an order of magnitude more MC 
steps to reach the same level of accuracy. The stopping criterion for the self-consistency 
loop was that two successive values of the static local susceptibility differed by less than 
0.5%. This is about twice the size of the statistical error. On the basis of these figures, we 
estimate that our final results for xt are accurate to within 1%. 

IV. RESULTS 

We have simulated the mean-field Kondo alloy model for fixed values of the Kondo cou- 
plings and several values of J for T > 0.05 T^ where T^ is the single-site Kondo temperature 
(see below). The first step in the calculation is the choice of the parameters a, Jf and tq that 
define the underlying single- impurity Kondo problem (cf. Eq. (P)). As the low-temperature 
properties of all antiferromagnetic Kondo models are described by the same fixed point, we 
are free to choose these parameters using criteria of numerical convenience. For the particu- 
lar case a = 1/2 and for all values of J±, the Kondo model is equivalent to a simple exactly 
solvable problem, the resonant model [^. We have made this choice as it provides us with 
means to test our numerical methods by comparison of the Monte Carlo results means with 
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the analytical solution. We have taken Tq ^ as the unit of energy and we have arbitrarily set 

Fig. |l] shows x(t) as a function of the scaled variable t//3 for J = and several tempera- 
tures. The correlation function in imaginary time is real and symmetric around r = /3/2 as a 
consequence of time-reversal invariance. Its minimum value steadily decreases with decreas- 
ing temperature. The expected behavior of the zero-temperature dynamic susceptibility in 
the long-time limit is xi'^) o^ ''""^ EH- At finite temperatures this expression generalizes to 



x{t) ~ (7r//?)^sin ^(vrr//?) ||2^. We have fitted our data for r ~ j3/2 with the expression 



\p J sm^ + sm^ 

where A is a T-dependent amplitude and the cutoff f is of the order of the inverse of the 
Kondo temperature to be defined below. The fits, shown by the solid lines in the figure, are 
in excellent agreement with the numerical data. 

The static spin susceptibility has been computed from the Monte Carlo results using the 
expression ||13| xt = Jo "TxI^)- The results thus obtained are shown in in Fig. |^. We also 
show for comparison the susceptibility of the resonant model, 

where (j){z) = (Plnr{z)/dz'^ and A is the width of the resonant level. The latter has been 



determined by fitting the data for tqT < 0.2 to Eq. ([T^) with the result A = 0.827 t^^. There 
is very good agreement between the theoretical expression and the Monte Carlo results in 
this temperature range. Deviations from the theoretical result are expected (and observed) 
at higher temperatures as Eq. (|18|) is only valid for Ttq ^ 1. Taking the T —>■ limit of 
Eq. ( [I8| ) we find the zero-temperature susceptibility xo = 2/(7rA) ^ 0.77 tq. Defining the 
Kondo temperature by xo = 1/(2 T^), we obtain T^tq ^ 0.65 . 

We have similarly computed the r-dependent susceptibility of the system for several 
values of J 7^ and T. The overall shape of the curves thus obtained is similar to that of 
those of Fig. ^ but the decay of the correlations becomes slower and slower as J increases. 

12 



This is shown in Fig. ^ where we show results obtained for several values of J at our lowest 
temperature, Ttq = 32^^. This slowing down of the spin dynamics, which is accompanied 
of an increase of the susceptibility, is a precursor effect of the phase transition that, as we 
shall see next, takes place for sufficiently large J . 

A necessary condition for the stability of the paramagnetic phase is that the inequality 
is Y{T) = 1 — Jxt > holds |[13|. Y{T) is plotted versus temperature in Fig. § for several 
values of J. The symbols are the Monte Carlo data. The dashed lines are fits of the results 
to a model that will be discussed below and that allows us to extrapolate the results down 
to T — > 0. We see that Y{T = 0) vanishes for J = 0.75 r^^ ~ 1.15 T^ which identifies it as 
the critical coupling. For J > J^ y{T) vanishes at a finite temperature, Tg. 

Before discussing in detail the temperature dependence of the uniform susceptibility in 
the vicinity of the critical coupling, we shall make a disgression in order to derive a simple 
model in terms of which the numerical data can be analyzed in a transparent way. We start 
by noticing that the frequency-dependent susceptibility may be related to the fiuctuations 
of the auxiliary field ^(t). Using Eqs. (^) and (|^) one may readily show that 

(l ^{i^n) 1^) = QiiuJn) [1 + x{^^n)Q{iuJn)] , (19) 

where the expectation value on the lefthand side of the equation is taken with respect to 
the probability distribution V[C,{t)] oc exp(— /3JF) with 

r/3 



:F[ar)]=m + ^T.W^-kTln ( Texp 







dTar)S,{T) 



(20) 



Here, pj^[0] = - In [2 cosh (/3r/2)] and ((...) )p = Tr [(. . .) exp (PTS^)] /Tr [exp {PVS^)]. 
Assuming for the moment that the transverse part of the effective field dominates over 
its fiuctuating longitudinal component, F ^ < ^'^(r) >^/^, the free-energy ( pOD may be 
expanded up to second order in ^ : 

J^l^r)] = ^[0] + ^E [ Q''M - X^'\u:n) ] U^n)? + . . . , (21) 

where the zeroth-order transverse susceptibility x^^\^n) = Tma;/ (F^ + cj^) and tti^ = 

13 



l/2tanh(/3r/2). Combining Eqs. (|2T|) and ([T9|) we derive the following expression for the 
frequency-dependent susceptibility : 



where -ft'(ct;„) is the Fourier transform of the first term on the righthand side of Eq. ([lOD. In 
the limit Ttq, |ti;„,|ro ^ 1 this is 

K{u;n) = {I -\uj„\ro + ...). (23) 

To 

Substituting this expansion in Eq. (0) we find : 



X{^n> = jj < — + Tft: + a\ujn\ - ^ 



u;2 



n I rri I ~ I I 

- + Ik + a|^n| 



' -jn (24) 



where a = na, Tk = T — o/tq and we have assumed T ^ F. From Eqs. (P^ ) and ( [T^ ) we 
can estimate (^^(r)) ^ (J^ + 2a /tq) at large F. We therefore expect Eq. (|2^ ) to be valid 
provided the condition 



Fro ^ max{ Jto, v2a}, (25) 



is satisfied. It is clear that Eq. ( P5| ) will not be fulfilled by the bare parameters, in general. 
We can, however, imagine writing down a set of renormalizat ion-group equations for the 
flow of the different couplings as the high energy cutoff is reduced. By analogy with the 
single-impurity case, we expect that in the paramagnetic phase the Kondo couplings will 
flow to the strong coupling fixed-point T = J^ —>■ oo. Therefore, we expect Eq. ( p4|) to 
become appropriate below some energy scale with renormalized values of the couplings. 

The derivation of the renormalization-group equations for the Kondo-alloy model is out- 
side the scope of this work |2^. We shall instead consider (^) as a phenomenological 
equation containing three renormalized parameters, T^, a, and F to be determined by a fit 
of the numerical results. 

The ratio Tk/J is determined by the static uniform susceptibility alone. Setting u;„ = 
in Eq. (|2^) we have : 
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■^2 



JXT = ^ - /^ - 1- (26) 

We thus see that the instabihty of the paramagnetic phase is signaled by the vanishing of 
the quantity under the square root in Eq. (pB]). We may therefore take A = T^/J"^ — 1 as a 
measure of the distance to the quantum critical point and rewrite the susceptibility in the 
form : 

JXt = Vl + A - VA. (27) 

Notice that, even if the assumptions made in the derivation of Eq. (^7D are not valid, 
the latter can still be regarded as a parameterization of the susceptibility in terms of a new 
quantity, A(T, J). The interest of this parameterization stems from the fact that the T- and 
J-dependence of A is very simple. We show in Fig. (|]) the numerical values of A obtained 
inserting the Monte Carlo results for the static susceptibility in Eq. (|27|). The dashed lines 
are fits to the simple functional form : 

A(T) = Ao + {T/Tof\ (28) 

where the parameters Aq and To are functions of J but not of T. The fits are very accurate 
over the entire temperature range of our simulations. The lowest curve, corresponding to 
our estimated value for the critical coupling, has been fitted with Aq = 0. Examination of 
the J dependence of Aq shows that, near Jc, Aq ^ a {1 — J/Jc), with a ^ 1. The parameter 
To has a finite limit. To ~ 0.27 T^ as J — ^ Jc- 

Eqs. ( p7D and (|28| ) imply that, in the neighborhood of the critical coupling, the uniform 
susceptibility has a non-Fermi-liquid T-dependence : 

'l-(S)''' for Af«T/To«l, 

JcXt ~ I (29) 

^l-V^-i^ter f- T/To«Af. 
The values of the exponent 7 found in the two regions defined above correspond to 
those obtained by other authors @-§l in the quantum critical (QC) and quantum disordered 
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I (QDI) regions in their analysis of metallic spin-glass models. A third region (QDII) is 
expected at lower temperatures where the normal Fermi-liquid T^ behavior is recovered. 
This second crossover is not visible in our data because, as we shall see, it occurs below the 
lowest temperature that we can reach in our quantum Monte Carlo simulations. 

The remaining two parameters in Eq. (^^ may be determined from an analysis of the 
full ojn- and T-dependence of the susceptibility. This is shown in Fig. ^ for two couplings, 
J = 0.65 Jc and J = Jc- The symbols are the quantum Monte Carlo data. The dashed 
lines are plots of Eq. (|2^) with F, a and Tk adjusted to fit the data. The quality of the fits 
is excellent for all the values J and T considered. It is remarkable that all our numerical 
results could be fitted with the same values of Ftq ~ 2.4 and a ~ 1.48. The full J- and T- 
dependence is therefore in the effective Kondo temperature, T^. The values of T^ obtained 
from the fits of the frequency dependence of the susceptibility are consistent with those 
determined above from the static susceptibility. 

At this point, we can discuss the connection between our results and those of Sengupta 
and Georges [^]. Phrased in our language, their approximation for xi^n) is obtained by 
ignoring the transverse field term in Eq. (|TT]), and solving the remaining Ising problem in 
the spherical approximation. The expression for the dynamic susceptibility that results from 
this procedure is : 



1 



Xi^n) ~ -^ j-^ + "kni - V [^ + Oi\uJn\]^ ~ "^ J ' ^^^^ 

where A is a Lagrange multiplier introduced to impose the spherical constraint, 

1 oo 

{S.{r)S^{r)) = - Y. XM = 1/4. (31) 

f-^ n=~oo 

As it stands, expression (^) diverges because Eq. ( pOD does not have the correct uj~^ high- 
frequency behavior. A high energy cutoff A must therefore be introduced in the simplified 
model. 

Comparison of Eqs. ( P^ and (^) shows that the two expressions become equivalent at 
low-frequencies provided we identify T^ with A and aT with A. The fact that the Monte 
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Carlo data could be fitted using Eq. ( p4D with T- and J-independent values of a and F 
justifies a 'posteriori the use of a constant cutoff in the simplified model. Once this is fixed, 
the parameter A can be determined from condition ([3l|). Since the numerical data do satisfy 
this normalization, it is not surprising that the determination of Tk from fits of the Monte 
Carlo data and that of A from enforcement of Eq. ( pi| ) result in the same temperature 
dependence. This suggests to use Eq. (|2^ ) in conjunction with the normalization condition 
to estimate the T-dependence of A in the temperature region for which we do not posses 
numerical data. We show in Fig. ^ the result of applying this procedure to the case of 
J = 0.8 Jc- We see that the expected crossover from a T^/^ law to normal T^ behavior 
occurs at a temperature T* ?a 0.06 T^. This is at the lower end of the temperature range 
that we can reach. The crossover temperature further diminishes as, J ^ Jc where it 
vanishes. This explains why normal behavior has not been seen in our simulations. 

The imaginary part of the magnetic response can now be determined by analytic con- 
tinuation of Eq. (p^. The general expression is complicated and not very illuminating. 
However, in the low frequency limit, cj <^ a F, and for J -^ J^, Eq. ( P^ can be cast in the 
scaling form : 

J.X"(-) = v^$ (^) , (32) 



where the universal scaling function $(x) 



A 1/2 



-1/2 



(l+a;2) +1 . (33) 



1 

This expression is equivalent to that found in references 

The low- frequency behavior of x"('^) as, J ^ Jc follows from the above equations : 

ff^ for cu < JcA, 

JcX\^) ^ I (34) 

(^)^^^ for JeA<cj<aF. 

The dissipative part of the susceptibility in the limit cj ^ is Fermi-liquid-like everywhere 
except at the quantum critical point. However, the characteristic spin-fluctuation frequency 
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ujsi oc JcA vanishes as J ^ Jc with cUsf oc (1 — J/Jc) for J ^ Jc and cjsf oc T^/^ at the 
critical couphng. The behavior of x"{^) ^^ Jc is non-Fermi-hquid-hke, x"(^) '^ \f^ which 
reflects the slow decay of the time-dependent spin-spin correlation function, < Szif)SziS^^ >~ 
t~^/^ that anticipates the appearence of long-range order in the system. The non-trivial T- 
dependence of the spin-fluctuation frequency in the vicinity of the quantum critical point 
that is responsible for the singular behavior of the susceptibility gives rise to anomalous 
powers in other thermodynamic and transport properties as well. In particular, Eqs. ( |32D 
and ( p3D imply that the temperature corrections to the specific heat and the resistivity 
behave, respectively, as 5C /T oc — vT and 5p oc T^/^, in the quantum critical region 0. 

The full cj-dependence of the absorptive part of the dynamic susceptibility is shown 
in Fig. ( H) for several temperatures at the critical coupling and several values of J at 
T = 0.05 T^. These curves have been computed by analytically continuing the fits of the 
imaginary-frequency Monte Carlo data. These curves are very similar in shape to those 
obtained in reference pO) for the single-impurity Kondo model and may be characterized 



by an effective Kondo temperature that decreases with the distance to the quantum critical 
point where it vanishes. Indeed, Eqs. (^) and (|33D imply that the effective Kondo scale uk-, 
defined as the half- width of the relaxation function x"(u;)/u;, is ujk ~ T\^ \/~K. 

V. CONCLUSIONS 

In this paper, we have studied numerically a Kondo lattice model with random exchange 
between localized spins. A mapping of this model to a self-consistent single-spin problem, 
exact in the limit of large lattice coordination, allowed us to obtain a complete numerical 
solution of the problem. The system has a quantum critical point between a normal metal 
and a spin-glass state. There is a region in the T — J plane near the quantum critical point 
where the characteristic spin-fluctuation energy varies as a non-trivial power of temperature 
(tUsf oc T^'^). This gives rise to non- Fermi-liquid behavior in thermodynamic and transport 
properties. At low enough temperature normal Fermi liquid behavior is recovered, except at 



the critical coupling. Our numerical results can be very well described over a large range of 
frequency and temperature by a simple model that we derive in the strong-coupling limit. 
This model is closely related to the M-component quantum-rotor and mean-field models that 
have been previously discussed in the literature @-§[. Some interesting questions remain 
open, notably, to what extent the assumption that the electronic bath remains unrenor- 
malized is a valid one. The enforcement of the self-consistency condition (^) poses some 
important technical difficulties which we hope to be able to overcome in future work. One of 
us (M. J. R.) acknowledges support of Fundacion Antorchas, CONICET (PID A^°4547/96), 
and ANPCYT (PMT-PICT1855). 
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FIGURES 
FIG. 1. Dynamic spin susceptibility in imaginary time for J = for (3/tq= 18, 20, 24, 28 and 

32, from top to bottom. The symbols are the Monte Carfo data. The sohd hnes are the fits referred 

to in the text. 

FIG. 2. Temperature dependence of the susceptibility for J = 0. The circles are the Monte 
Carlo results. The dashed line represents the susceptibility of the resonant model. 

FIG. 3. The J-dependence of the dynamic spin susceptibility in imaginary time at fixed tem- 
perature, T To = 1/32. The values of Jtq are 0, 0.6, 0.65, 0.7, 0.75, from bottom to top. The symbols 
are the Monte Carlo data. The lines are guides for the eye. 

FIG. 4. The J- and T-dependence of Y[J, T) = 1 — JxT- The paramagnetic phase is stable for 
Y{J,T) >0. 

FIG. 5. The J- and T-dependence of A. The symbols are the Monte Carlo data. The dashed 
lines are the fits to the expression in Eq. (pq). 

FIG. 6. The to- and T-dependence of the dynamic susceptibility for two values of the coupling 



J. The symbols are the Monte Carlo data. The lines are fits to the expression in Eq. (24). 



FIG. 7. The crossover between anomalous and Fermi-liquid behavior obtained from Eq. (24) 
and condition (|3l|). 



FIG. 8. The relaxation function, x"(w)/w. The curves have been obtained by analytic contin- 
uation of the fits of the imaginary-time Monte Carlo data, (a) T/Tk = 0.05 and Jtq = 0.55, 0.6, 
0.65, 0.7 and 0.75, from top to bottom, (b) J = Jc and (3/tq = 32, 28, 24, 20, and 16, from top to 
bottom. 
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